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A closed-form solution to the energy-based stochastic Schrodinger equation with 
q ; a time-dependent coupling is obtained. The solution is algebraic in character, and is 

' expressed directly in terms of independent random data. The data consist of (i) a ran- 

dom variable H which has the distribution ¥(H = Ej) = 7Tj, where 7Tj is the transition 
probability KV'ol^j)! 2 from the initial state |V>o) to the Liiders state \<pi) with energy 
Ei\ and (ii) an independent P-Brownian motion, where P is the physical probability 
measure associated with the dynamics of the reduction process. When the coupling 
(T) ■ is time- independent, it is known that state reduction occurs asymptotically — that is 

to say, over an infinite time horizon. In the case of a time-dependent coupling, we 
■^j- . show that if the magnitude of the coupling decreases sufficiently rapidly, then the 

energy variance will be reduced under the dynamics, but the state need not reach an 
energy eigenstate. This situation corresponds to the case of a "partial" or "incom- 
plete" measurement of the energy. We also construct an example of a model where 
the opposite situation prevails, in which complete state reduction is achieved after 
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the passage of a finite period of time. 



PACS numbers: 03.65.Ta, 02.50.Ey, 02.50.Ga, 02,50,CN 

I. INTRODUCTION 



^ ■ This paper is concerned with the problem of obtaining closed-form solutions to the energy- 

based stochastic extension of the Schrodinger equation in the case of a time- dependent cou- 
pling parameter. In this situation the dynamical equation of the wave function is assumed 
to satisfy the following stochastic differential equation: 

d|V t ) = -iH\ipt)dt - \a 2 t {H - H t ) 2 \^ t )dt + \a t [H - H t )\ijj t )dW t . (1) 

Here {\ipt)}o<t<oo denotes the state- vector process, which is defined on a probability space 
(fi,jF, P) with filtration {J r t}o<t<oo, with respect to which {W t }o<t<oo is a standard one- 
dimensional Brownian motion. The operator H is the Hamiltonian of the system, and 

H t = (^ t \H\^ t ) (2) 

is the expectation value of H in the state \ipt)- The time-dependent coupling parameter 
{ct}o<i<oo, which has the units [energy] -1 [time] ~ 1//2 , and is assumed to be a given positive 
function, determines the rate at which state-vector reduction occurs. For simplicity we 
consider the case of a pure initial state iV'o); the generalisation to the case of a mixed initial 
state is straightforward. For convenience we assume that the initial state has norm unity; a 
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straightforward exercise making use of the Ito calculus then shows that (ip t \ipt) = 1 for all 
t > 0. Likewise for simplicity we assume that the Hamiltonian H has a discrete spectrum, 
and that the Hilbert space is of finite dimension. We assume that the possible values of 
the energy are given by {Ei\i = \^...,N- The transition probability from the initial state j^o) 
to a state of energy Ei is Hi = (^olfiil^o), where IT denotes the projection operator onto 
the Hilbert subspace of states with energy Ei. Equivalently, we can write 7Tj = | (</>i|^o) | 2 , 

— 1/2 A 

where \<pi) = n i Hi\ipo) is the so-called Liiders state associated with the given initial state 
|"0o) and the energy E^. We note that (6Ad>i) = 1. According to the von Neumann- Liiders 
state- vector reduction hypothesis 0,EEEo|, if the initial state of the quantum system is |^o) 
and a measurement of the energy is made with the result E iy then there is a discontinuous 
transformation of the state of the system, and the new state is given by \<fii). In particular, 
in the case of a degenerate spectrum a specific state is selected in this way among all those 
with the given eigenvalue. 

In the stochastic framework the discontinuous von Neumann-Luders reduction process is 
replaced by a continuous reduction process modelled by equation Depending on the 
details of the physical setup, the reduction process modelled by (JTJ) can be regarded as taking 
place either (a) as a consequence of a measurement having been made, or (b) as a result of 
interaction of the system with its environment, or (c) spontaneously. For instance, we can 
view (0) as a phenomenological "reduced-form" model for the dynamics of a system when 
an energy measurement is made. The time-dependent coupling in that case represents the 
exogenous intervention of the measurement apparatus in the dynamics of the system. 

The mathematical and physical properties of the energy-based stochastic extension of the 
Schrodinger equation have been studied extensively in the literature [H 0, H 0, El EH EH 
121 . 13L 22L 24l 32l . 33l | . The energy-based stochastic Schrodinger equation can be regarded as a 



special case of a more general class of stochastic models for the dynamics of the wave function 
that use nonlinear stochastic differential equations of the form (^Q), but typically involving 
a number of quantum operators. Such models have been introduced with a variety of aims, 
and in the general situation the operators need not commute amongst themselves or with the 
Hamiltonian, and need not be Hermitian. We mention, e.g., B HE EE Ei H El H 0, 
and works cited therein, as references to the substantial body of publications in this area. 

The significance of the energy-based reduction model, in contrast with the more general 
situation, is that energy is conserved in expectation. More specifically, it is a property of the 
dynamical equation ((TJ) that ti{p t H) is constant, where the time-dependent density matrix p t 
is defined by p t = K p [\ilj t )(ip t \}. Here E p [— ] denotes expectation with respect to the physical 
probability measure P. Conservation of the energy follows from the fact that 

^ = -i[H, p t ) + \a\ [Hp t H - \H*p t - \p t fr] . (3) 

Thus, the energy-based stochastic Schrodinger equation can be regarded as appropriate to 
the description of the dynamics of an isolated system, or any system for which on average 
there is no net exchange of energy with the environment. 

When the coupling parameter in equation (JJJ is a constant, it ispossible to obtain a 
closed- form solution to this stochastic differential equation U, E3, 13 • The solution 



to Q in the case of a constant coupling parameter can at each time t be expressed as a 
function of a state variable £t, the value of which is determined by the specification of a pair of 
independent random data corresponding roughly to the idea of a split between "signal" and 
"noise". More specifically, the state-variable process {^}o<t<oo is of the form £ t = o~Ht + B t , 
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where if is a random variable that takes the value Ei (i = 1, 2, . . . , N) with probability 7Tj, 
and {B t } < t<oo is an independent Brownian motion. The purpose of this paper is to show 
how this construction can be generalised to the case of a generic time-dependent coupling. 

The paper is organised as follows. In Section |H] we introduce an ansatz that leads to the 
solution of Some comments are made on the interpretation of the ansatz and its relation 
to similar constructions in the theory of filtering. In Section lTTTl we obtain general expressions 
for the conditional probability 7r# = E[l {H=E i }\J r t\ and the conditional energy expectation 
H t = E[iJ|jF t ] in terms of the process {£t}o<t<oo- Whereas in the time-independent case the 
random variable H t can be expressed, for each value of t, as a function of the state variable 
£t, in the time-dependent case H t is a functional of the trajectory {^ s }o<s<t- In Section ITVl 
we show that although {£ t } is in general non-Markovian, the corresponding energy process 
{H t }, which is a functional of does have the Markov property. 

In Section we show how the analysis of Section IIHI leads to an expression for the 
state- vector process that solves (Q). In Section IVTI we examine the corresponding "inverse" 
problem. We show how, given any random trajectory \ip t ) that solves (Q), it is possible 
to construct explicitly the independent random data H and {-£><} associated with it. We 
also remark on the relation between our solution technique and a well-known linearisation 
method often used in analysing the properties of Then in Section [VIII we show that 
the state vector collapses to the specified eigenstate, provided a condition on the coupling 
is satisfied. 

More generally, since we consider an essentially arbitrary time-dependent coupling, one 
can envisage circumstances in which the coupling ceases before state reduction is complete. 
This kind of situation can be regarded as representing an approximate measurement of 
the energy, in which partial information is gained but no definite outcome is obtained. In 
Section IVIIII we examine this case in some detail, and derive upper and lower bounds for 
the asymptotic value of the energy variance as t goes to infinity. 

Another interesting situation arises when the magnitude of the coupling increases suf- 
ficiently rapidly to ensure that state reduction is completed after the passage of a finite 
amount of time. In Section IIXI a special example of such a model is constructed. This ex- 
ample turns out to have a direct relationship to the finite-time collapse model introduced in 
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The model presented in [11[ is based on a constant coupling parameter and a Brownian 
bridge noise. Here we present a alternative model for finite-time collapse, for which the noise 
is a standard Brownian motion and the coupling is time- dependent. We demonstrate that 
the two models are physically equivalent. 

While the probabilistic method presented in Section ITTT1 is effective in obtaining the so- 
lution to the problem under consideration here, there are circumstances in which other 
methods are useful as well. We thus outline two further approaches for obtaining the solu- 
tion. In Appendix El we introduce a method that involves a discretisation and a continuum 
limit, and proceeds in a manner similar to the analysis entailed in the evaluation of Feynman 
integrals. The method is computationally intense, but is satisfying because it allows one to 
work directly with the quantities under investigation. In Appendix [B] we consider another 
method that is similar to the path integral approach, except that we use a decomposition 
of the state-variable trajectory into increments, and we regard the time-dependent coupling 
as moderating the noise rather than the signal. 
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II. THE QUANTUM INFORMATION PROCESS 



The ansatz that we use to solve (JTJ) is based on the specification of a state- variable process 
{£t} which, for reasons discussed shortly, will be called the "quantum information process", 
and is of the following form: 



Here H denotes a random variable on the given probability space (fl,F, P), taking the 
possible values {-E'i}i=i,2,...,v with the probabilities {iri}i=i,2,...,N, an d {B t }o<t<oo is a standard 
Brownian motion, independent of H. We do not assume that {B t } is adapted to the filtration 
{J-'t} introduced earlier. On the contrary, we shall see later that {Ft} is generated by {&}. 
The various terms appearing in (£Q) can be given an interpretation in the language of filtering 
theory. This "signal and noise" interpretation, although not essential to the use of the 
ansatz to solve (|TJ), is nonetheless physically very suggestive, and as a consequence helps to 
motivate the form that the solution takes. Indeed, the methodology of filtering theory has 
been already shown to be effective in deriving solutions to the energy-based stochastic 
Schrodinger equation, and in what follows we take this line of investigation further. 

The random variable H, according to this interpretation, represents the unknown terminal 
value of the energy of the quantum system whose time evolution is described by equation (JTJ). 
The term H J Q * a s ds in (j3J) should be thought of as the "signal" component of the quantum 
information process. As time passes, the magnitude of the signal component increases, but 
the true value of H remains obscured by the presence of a noise process {B t }o< t <oo- The 
"accessible" information concerning the value of H is represented by the process {^}o<t<oo, 
which consists of both signal and noise. Given the history {£ s }o<s<t over a finite time interval 
[0, t], it is not generally possible to disentangle the true value of H from the noise. In the 
context of filtering theory, the task in such a setup is to determine the best estimate of H, 
given the information of the trajectory {£ s }o<s<t from time zero to time t. It is a remarkable 
feature of the stochastic Schrodinger equation that the expectation of the Hamiltonian turns 
out to be given by such an estimate. 

In the first part of the paper we shall examine the case for which the state- vector trajectory 
{iV't)} satisfies the dynamical equation (JTJ) for all t e [0, oo). In order for the trajectory 
{|^t)}o<t<oo to be well-defined, the coupling function {a t }o<t<oo must be such that 



Additionally, as will be established in Section I VI II to ensure a complete reduction of the 
state vector, the coupling function must be chosen such that 



The purpose of this condition is to ensure that the coupling remains reasonably "strong" 
for all time, and does not attenuate too much. 

There are circumstances in which (JBJ) is satisfied but © is not. In such situations the 
wave function need not fully collapse to an eigenstate, even though the energy variance will 
be reduced. This case is examined in Section IVIIII It will be assumed throughout Sections 




(4) 




< t < oo. 



(5) 




(6) 
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IIimVHll that {£ t } and {cr t } are defined for all t in the range < t < oo, and that (0) holds. 
In Section HXl we drop the assumption of an infinite collapse time, and consider the case for 
which the integral of {a t } diverges after a finite passage of time. 

The estimation problem posed by an ansatz of the form (J3J) is well established in the lit- 
erature of nonlinear filtering [HQ. The relevance of (gj) to the dynamics of the quantum 
state {l^t}} satisfying the stochastic Schrodinger equation ((TJ), on the other hand, is not 
obvious. As we shall demonstrate, the information generated by {£ s }o<s<t is equivalent to 
the information generated by the evolution {\ip s )}o<s<t of the quantum state itself. We for- 
malise this notion by observing that {^}o<*<oo and {\?pt)}o<t<oo generate the same filtration 
{J-'t}o<t<oo- As a consequence, the energy process {H t } determined by the quantum expec- 
tation (J2J) of the Hamiltonian operator turns out to be indistinguishable from the process 
generated as t varies by the mathematical expectation of the random variable H, conditional 
on the specification of the trajectory {£ s }o< s <t- 

The best estimate for H, in the sense of least quadratic error, given the history of the 
information process up to time t, is known (see, e.g., e lGj) to be the conditional expectation 

F t = E[#|{£ s } Q < s < t ]. (7) 

We have used the same notation {H t } for the processes defined in (0) and in (J7J) because 
these processes will be shown to be the same. When {a t } is constant, equation (J7J) can be 
simplified to the form H+ = E [-ff|£t]- In this case, {£ t } is Markovian: this is the situation 
considered in 

0, Q H El • However, if {a t } is not constant, then in general the trajectory 



{£ s }o<s<i must be taken into account to determine the conditional expectation ((7|). The non- 
Markovian property of {^} can be seen intuitively as follows. Writing (0J) in differential form, 
we have 

d& = (T t Hdt + dB t , (8) 

which makes it evident that {at} determines the strength of the signal, that is to say, the 
rate at which the true value of H is revealed. If {at} is constant, then sampling from 
at any small time period in the interval [0, t] is as good as any other. This is, in essence, the 
property of when {a t } is constant, that makes it Markovian. If {a t } is not constant, 
then there is a temporal bias in the sampling from and observations from different 

periods cannot be treated on an equal footing. 



III. CONDITIONAL PROBABILITY PROCESS 

Our goal in this section is to work out an explicit expression for the conditional probability 
process {ttu} defined by 

7T ft = E[l {lf=Jg . } |{&}o<,<t]- (9) 

Here denotes the indicator function which takes the value unity if A is true, and zero 
if A is false. Once we obtain 7i it , then the conditional expectation (J7J) can be obtained by 
the relation H t = ir it Ei. 

To determine {71"^} we use a change-of-measure technique (see, e.g., jl^, 17, 26[), proceed- 



ing as follows. Let (Q, J 7 , P) be a probability space on which a standard Brownian motion 
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{B t } is defined, and let if be a random variable on (f2,jF, P) that is independent of {B t }. 
We fix a time interval [0, u], and let {£t}o<t<« be given by 

& = H I a s ds + B t , (10) 
Jo 

where {a t } is deterministic and satisfies (jHJ)- Next, we define a process {A t } < t < u over [0, u) 
by the expression 

At = exp (h J a s d£s - \H 2 jf cr s 2 d^ , (11) 
or equivalently, by virtue of (fTUJl . 

A" 1 = exp j£ - <^ ds l • ( 12 ) 

The idea is to use {At} to make a change of probability measure. The new probability 
measure will be defined on the space (fl,Q u ), where Q u C T is the a-subalgebra of events 
determined by the specification of the trajectory {B t } < t < u over the given time horizon, 
together with H. 

We recall that the points of Q represent the possible outcomes of chance, and the elements 
of T are subsets of Q with the property that for each such subset A G T the measure P 
assigns a probability F(A) to the event that uj G A. The elements of Q u consist of those 
elements A G T with the property that knowledge of the value of H and the trajectory 
{B t }o<t<u is sufficient to determine whether u G A. For any element A G Q u the value of 
the indicator function 1{u}<=a} is determined by the specification of H and {B t } < t < u . The 
new measure Q on (fl, Q u ) is then given as follows. For any set A G Q u we define 

Q(A)=E p [A- 1 l {weA} ]. (13) 

This relation is usually abbreviated by writing dQ = A~ x dP. Since E p [A" 1 ] = 1 by virtue 
of elementary properties of the stochastic exponential (|12|). it follows that Q(Q) = 1. 

Given the setup described above, we have the following facts: (i) on the probability space 
(ft, Q u , Q), the process {£t}o<t<« defined by (fTUj) is a Brownian motion, and is independent 
of H; (ii) the random variable H has the same probability law with respect to Q as it does 
with respect to P; (iii) for all t G [0,u], the conditional expectation f t = K ¥ [f '(if )|{£ s } < s <t] 
of a function of the random variable H can be expressed in the form 

_ EQ[/(ff)A f |{6,}o<,<t] (u] 
Jt EQ[A t |{e s } < s <t] ' 1 ] 

To work out an expression for T\ it we start with the definition substitute (|TT|) into 
(JHJ) and set /(if) = 1 {_h"=_b i; } to obtain 

7Ti exp /„* a s d£s - \Ef jj (rfds ) 

TTit = ^ i '—. (15) 

Ei ^ ex P ( E * Jo a sd^ s - \E} J Q a 2 s ds 
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Similarly, by setting f(H) = H we obtain 

Ei *i E i ex P /o CT * d ^ ~ \ E i Jo a s ds 



Ei TTi exp ( /„* <7 s d£ s - |£ 2 /„* a 2 As 



(16) 



The result appears at first glance to depend on the choice of time horizon u, since {ttu} and 
{H t } are only defined for t G [0, «]; but it is straightforward to see that (fTK|) and (fTHjl remain 
valid for all i G [0, oo). 

The arguments establishing the validity of (i), (ii), and (iii) above can be sketched briefly 
as follows. First we note that the relation between expectation under the P-measure and 
expectation under the Q-measure is given by 



E Q [X] = E V [A- 1 X]. 
To see that {£i}o<t<oo is a Brownian motion on (fl,J-, 



it suffices to show that 



±(x 2 s+y 2 t+2xys) 



(17) 



(18) 



for s < t and for all x, y. In particular, if {£ t } possesses this bi-characteristic function, then it 
follows at once that {£ t } is Q-Gaussian, and that Cov[£ s ,£ t ] = s for s <t. These properties, 
together with the fact that is continuous, are sufficient to characterise it as a Brownian 
motion under Q. The verification of (fT%|) follows by a calculation that makes use of fllOL 
f!12jl . and (fT7j) . the P-independence of H and {B t }, and basic properties of {B t } under P. 
The Q- independence of H and {^} then follows by a similar calculation that establishes 
that 



71; C 



(19) 



To show that the probability law of H is the same under P and 
P-independence of H and {B t } to observe that 



we make use of the 



\H = E l ) = E p [A~ 1 l{ H=Ei }] 



E 



E 



l{H= £l } exp 
l{//=B l} exp 



H / aAB„ 



EF [l {H=Ei} ] E* 
P (H = . 



exp 



(7 



4B X 



\H 2 



1 p 2 



a 2 ds 



cr ds 



crAB s 



1 7712 

2^ 



<r 2 ds 



(20) 



As for (iii), we remark first that relation ()14J1 is a special case of a more general result 
referred to as the Kallianpur-Striebel formula [27j. A derivation of (pPfj) can be sketched as 
follows. We reverse the construction above and start with a probability space (Q, Q) on 
which {£t}o<t<oo is a standard Brownian motion, and H an independent random variable 
taking the values {-Ej}i=i i 2,...,v with the probabilities {^i}i=i^,...,N- We fix a time interval 
[0, u] and let Q u denote the cr-subalgebra of events generated by {£,t}o<t<u and H. Assuming 
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that {a t } satisfies (JSJ), we define the process {A t } < t < u by (fTTj). as before, and we introduce 
the measure P by setting dP = A u dQ, or equivalently, ¥(A) = E Q [A u l {ujeA} ] • Then the 
process {B t } < t < u defined by B t = £ t — H J*cr s ds is a Brownian motion with respect to P, 
and is P-independent of H. It is worthwhile emphasising that Tt = a {{^s}o<s<t)i whereas 
Q t = a (H, {t, s } <s<t) = cr(H,{B s } < s <t)- Thus T t C Gt- One can think of {F t } as the 
filtration generated by the dynamics of the state vector process {\ipt)}o<t<oo', whereas in the 
larger "nonphysical" filtration {Q t } the value of H is already "known" at time 0. It follows 
by the conditional form of the change of measure relation that 

rfr f /mirf i i _ E<Q [f(H)A u \{£ s } < s < t ] / 91 s 

" EQ[A u |{6}o<-<d • (21) 
Finally, we use the fact that if and {£ s }o<s<m are Q-independent to deduce that 

[f(H)A u \{i s } < s < t ] = E Q [f(H)A t \{^ s } < s < t } (22) 
for any choice of f(H), and as a consequence we deduce (fTlj). 



IV. ON THE MARKOVIAN NATURE OF THE ENERGY PROCESS 

Before we verify that our expression (fTKj) for the conditional expectation of H agrees 
with the energy process (j2J), it will be useful to show that {H t }, as defined by (|TH|) . has the 
Markov property. We note in particular that by virtue of (fT7)|) we have 



J2i KiE t exp [Eith - \Ef f* a 2 s ds 

H t = ^ ^, (23) 

E< n i ex P {EiVt - \E1 j <r 2 s ds) 



where the process {r] t }o<t<oo is defined by 



Vt = / °sd£ s . (24) 



o 

To show that {H t } is a Markov process it will suffice if we can show: (i) that {r] t } is a 
Markov process; and (ii) that rjt can be expressed as a function of H t , i.e. that the relation 
between rj t and H t is invertible. In particular, if the relation between rj t and H is invertible, 
then the filtration generated by {77^} is the same as the filtration generated by {Ht}. 

To verify that {i] t } is Markovian we must demonstrate that, for all T > t and for all 
x G M, we have 

< x\{T] a } < B < t ) = F^T < x\ Vt ). (25) 

Alternatively, it suffices to show that 

F(rj t < x\r) s ,r) Sl ,7) S2 , . . .,rj Sk ) =F(rj t < x\r) s ) (26) 

for any collection of times t, s, si, S2, ■ ■ ■ , Sk such that t > s > s± > S2 > • • • > Sk > 0. To 
check that these conditions are satisfied we proceed as follows. First we note that by virtue 
of (J24|) we have 

r h = H [ a 2 s ds + / <T 8 dB a , (27) 
Jo Jo 
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and therefore 

Vu Vv __ Jo^sdBs £ a s dB s 



Io a s ds Io a s ds Io a s ds Io a s ds 

for all u > v > 0. We shall establish that the process {f u }u>o defined by 



(28) 



f " a s dB s , 
f u = Jo r ' 29 

Jo a s ds 

appearing in (|28|) has independent increments. In particular, since {f u } is a Gaussian process 
it suffices to show that <fb — fa and fd — <p c are independent for all d > c > b > a > 0. To 
see this, we note that since <pb — f a and (fid ~ fc are Gaussian random variables, for their 
independence it is sufficient to verify that the covariance E[(<p& — (p a ){{p d — ip c )] vanishes. 
But this follows after a short calculation making use of the Wiener-Ito isometry 



E 



aAB s / (TAB. 



u/\v 

2 



a 2 s ds } (30) 



where u A v — min(tt, v). Similarly, one can verify that ip c is independent of (fib — fa for all 
c > b > a > 0. This follows from the fact that the increment v?oo — fc is independent of 
fb ~ fa- Next we observe that as a consequence of the definitions of {rjt} and {ft} we have 

¥(r] t < x\r] s ,r} Sl1 r) S2 ,...,r) Sk ) = P (r} t < x \r] s , <p a - f Sl1 f Sl - f S2 , . . . , (p,^ - f Sk ) ■ (31) 

However, since Tj t and rj s are independent of ip s — <p si , (p Sl — <p S2 , . . ., f Sk „ 1 — fs k , the desired 
result follows. 

To show that H t is invertible as a function of rjt it will suffice to show that for each fixed 
t the function 



J2i nEi exp (Eitj - \E1 f a 2 s ds 

H(v,t) = f ^ (32) 

Ei ex P ( E iV - \ E 1 fo °s ds ) 



is monotonic in ij. But this can be seen immediately, since 

dH(r], t) E< - H{rj, t)f exp (Eft - ±£? f* a 2 s ds 



drj 



Ei ^ ex P - 1^ 2 /o a2ds 



(33) 



which is positive for all values of 77. As a consequence we see that the process {H t } as 
defined by has the Markov property. In particular, we have the following equalities: 
P {H T < x\{H s } < s < t ) = P {H T < x\H t ) =F{H T < x\r) t ) =F{H T < x\F t ). 

It is interesting to note that the conditional probability process {7^} given by ([T5|) can 
be expressed in the form 

7r is exp (Et f* a u d£u - \E\ f o*du) 

Kit = K (34) 

J2i K is exp (Ei j s a u di u - \Ef j s a 2 u du ^ 
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The interpretation of this relation is that the energy-based reduction models exhibit a fun- 
damental dynamic consistency property. In other words, if at some intermediate time s, 
where < s < t, we take note of the a posteriori conditional probability ir is , which is based 
on all information available up to time s, then we see that the resulting "new" model for 
the collapse of the wave function, given by ()34j) . is of exactly the same form as the original 
model, with 7ii S playing the role of the new a priori probability This means that the choice 
of the initial time t — has no preferential status in the theory. Indeed, the dynamical 
consistency of the energy-based reduction theory shows that the objections raised by Pearle 
in this connection are groundless. It is perfectly consistent to regard the collapse process 
as having already started at some earlier time than "the present" . 



V. INNOVATION PROCESS AND SOLUTION 

In Section 11111 we calculated the expectation of the random variable H conditional on 
the specification of {£ s }o<s<t, and we claimed that the result gives the energy expectation 
process (J2J). The aim of this section is to verify this claim, by showing how {£ t } is related 
to the Brownian motion {Wt} of equation ([!}. We begin by analysing the dynamics of the 
energy process (j!6|) . A direct application of Ito's rule shows that 

dH t = -a 2 t H t V t dt + a t Vid£t, (35) 

where 

V t = E[(H-H t ) 2 \{£ s } < s < t ] (36) 
is the conditional variance of H . The next step is to define a random process {Wt} by 

W t = tt- [ a s H s ds. (37) 
Jo 

It follows then from (J35j) that the dynamical equation for {H t } is given by: 

dH t = <r t V t dW t . (38) 

Equation (|38|) can be given a simple heuristic interpretation if we write it in the form 

H t+dt - E [#t+d*|{6}o<.<t] = otVtdWt. (39) 

Since at t the variance Vt is known, as is also the conditional expectation ¥,[H t+ dt\{C,s}o<s<t], 
we see that dWt embodies the "new information" entering the system between t and t + dt. 
It is for this reason that {Wt} is called an innovation process. 

We claim that {Wt} is an {jF t }-Brownian motion. Here we recall that {Ft} denotes the 
filtration generated by the process {&}. Thus, conditioning with respect to the cr-algebra T t 
means conditioning with respect to the trajectory {£ s }o<s<t- To proceed we need now a more 
precise definition of Brownian motion. A process {W t } < t <oo on a probability space (Q, J 7 , P) 
with filtration {J r t}o<t<oo is said to be a standard Brownian motion if it satisfies the following 
properties: (i) Wo = almost surely; (ii) {W t } is {jF t }-adapted; and (iii) for all < s < t 
the increment Wt — W s is normally distributed with mean zero and variance t — s, and is 
independent of T s - In the present context we shall use the so-called Levy's characterisation 
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of Brownian motion, which states that if {W t } is a martingale, and if (dW t ) 2 = dt, then 
{W t } is a Brownian motion. 

Let us consider the martingale condition. Writing E t [— ] = E[— for conditional expec- 
tation with respect to T t , we shall establish that E t [Wr] = W t for t <T. We find 



a.HAs 



E t [W T ] = Et [fr] - E t 

= H t f a s ds + Et [St] - / <r s E t [H s ] ds, 
Jo Jo 



(40) 



where we have substituted and we have interchanged the order of integration and expec- 
tation by use of the Fubini theorem. Next, we note that 



I a s E t [H s ] ds = [ a s E t [H a ] ds + f a s E t [H s ] ds 
Jo Jo Jt 



cr s H s ds + H t a s ds. 



(41) 



Here we have used the fact that H t = E t [H] satisfies the martingale condition E t [i? s ] = H t 
for s > t. Hence substituting (j4T|) into (j4T)j) we obtain 



Et [W T ] =H t f cr s ds + Et [B T ] - [ a s H s ds. 
Jo Jo 

Finally, from the tower property of conditional expectation we have 

E t [B T ] = E t [l [B T \{B s } < s <t, H]j=E t [B t ] . 
Inserting this relation into (|42|) we obtain 



(42) 



(43) 



E t [W 7 



E t 
6- 



H / a s ds + B t 



cTeHAs 



a,HAs = W t 



(44) 



where we have used the relation E t [£ t ] = £ t . This establishes that {W t } is an {jFf}-martingale. 
On the other hand, because 



dW t = {H- H t )a t dt + dB 



t- 



(45) 



it follows that (dH^) 2 = dt. Taking this together with the fact that {W t } is an {T t }- 
martingale, we conclude by Levy's criterion that {Wt} is an {jF t }-Brownian motion. 

We are now closer to establishing the relation between (J2J and (|16j). To this end we 
consider the conditional probability (|T5j) that H takes the value J5j. Taking the stochastic 
differential of (|T3jl and substituting ([3*7jl into the result, we find, after some rearrangement, 
that {ttu} satisfies 



dn it = a t (Ei - Ht)n it dW t . 



(46) 
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With another application of the Ito formula, we thus deduce that 

d4 /2 = -\a 2 {E t - H t )\l/ 2 dt + \a t {Ei - H t )^dW t . (47) 

Finally, if we let denote the normalised Liiders state |H,|3(j associated with the eigenvalue 
Ei, and define {|^)} according to 

|^) = ^e-^4 /2 |0 4 ), (48) 

then it follows at once from (J47|) that {|'0t)} satisfies the time-dependent energy-based 
stochastic Schrodinger equation ((TJ). 

In summary, if we define {£t} and {t] t } in terms of the independent random data H and 
{B t } according to (gj) and ||23jl, and if we define {H+} and {[*/>*) } bv (fToTl and and {W t } 
by (jHZj)? then solves (JTJ) for the given Hamiltonian H and initial condition j^o)- 



VI. DIRECT DERIVATION OF INDEPENDENT RANDOM DATA 

The way in which we have solved equation is by introducing the concept of a quantum 
information process (gj) specified in terms of a random variable H and an independent 
Brownian motion {Bt}. It is possible, however, to deduce the existence of these random 
data directly from (0). In this section we shall illustrate this reverse construction. 

We begin by remarking that the energy-based stochastic Schrodinger equation can be 
cast into integral form, incorporating the initial condition jV'o)) as follows: 

\ipt) = exp (-iHt -\^a 2 (H- H s ) 2 ds + \j^a s {H- H s )dw)j |^o)- (49) 
After some simple rearrangement we then deduce that 

exp (-iHt + \H /' a s (dW s + a s H s ds) - \H 2 £ a 2 s ds) |^o) 

lV»t) = " 7 • (5°) 

exp ( \ /* a s H s (dW s + a s H s ds) - \ a^ds 



Given the {jF t }-adapted Brownian motion {W t } and the energy expectation process {H t } = 
{ip t \H\ip t ) we now define a process {^} by writing 

it = W t + f a s H s ds. (51) 
Jo 

It follows then that \ifj t ) can be written in the form 

exp (-iHt + \H f* a s dZ t - \H 2 jj a 2 s ds) |Vo) 
\A) = 7 ^ • (52) 

With these ingredients at hand we now claim the following: The random variables H = 
lim^oo H t and B t — £ t — H j Q a s ds are independent. Furthermore, the process {B t } thus 
defined is a standard F-Brownian motion. 
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The existence of the random variable H is ensured by the martingale convergence theorem. 
The fact that H has the distribution IP (if = E{) = 7Tj then follows as a consequence of known 
properties of the stochastic equation To show that H and B t are independent (for any 
value of t) it suffices to show that 



(53) 



for any x, y. The proof proceeds as follows. First, by use of the tower property of conditional 
expectation we have 



= E p 



Ab- H Jo °sds)+yH 



e! 



(54) 



where E f denotes conditional expectation with respect to Here we have used the fact 
that £ f is ^-measurable. For the inner expectation we can write 



E 



*(y- x fo °sds)H 



(y-x J Q a s As)Ei 



(55) 



where n it is defined as the conditional probability 

7i it = F{H = Ei\F t ). 



(56) 



To obtain an expression for {iru} we recall jj, |6J] that since the projection operator IT 
commutes with the Hamiltonian, the bounded process {(■?/>t|rLjV ; i)} is a martingale. It 
follows that 



(i/>t\fLi\iJ>t) = < (^ooin^oo) 



(57) 



On the other hand, by known properties of the reduction process (JTJ) we have ( < 0oo|nj|'0oo) 
1^=^}, and hence 



ME^) = Ef [l {H=Ei} ] 



(58) 



We are therefore able to deduce that P(i? = Ei\J-" t ) = Ha = (i/j t \tli\ip t ) . A short calculation 
making use of (J52j) and properties of the projection operator then shows that 



Tii exp [Ei jl a s d£ t - \Ef f* a 2 s ds 
exp ( Jl a s H s &i t - \ J* o 2 s H 2 s &s 



(59) 



In fact, with a little work one can also show that ttu = \(4>i\ijjt)\ 2 ■ In other words, 7r# is 
given by the usual formula for the quantum-mechanical transition probability from the state 
\ipt) to the Liiders state \4>i). Of course, in standard quantum mechanics it is an assumption 
that Ku, when defined in this way, has the interpretation of a transition probability. But in 
the stochastic theory, we deduce this property. 
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Returning to our calculation of the inner conditional expectation in formula ()54|) we see 
as a consequence of that 



It follows that 



V . 7T- e^ -a! Jo CT = ds ) B « e £! /o ^ d ?*-| E ? /o ^ 

e / t tf B ff s d£ t -±J t *ftffda ' 



e x & V- 7T- e^ _a: ^o °sds)-Ei Jo o- s dft-5-B? /„ o-^ds 
g/o'ff.fl-.dft-^/o'^Hfda 



(60) 



(61) 



Next we observe that since d£ 4 = dW t + a t H t dt, the expectation appearing above can be 
written in the form 



E 



e *$t e Ei Jo <r s d^-|B? J v 2 ds 



.(62) 



However, the expression 

$ t = exp 



(63) 



is the change-of-measure density over the interval [0, t] needed to make {£ s }o<s<t a standard 
Brownian motion. Writing Q for the resulting new measure, by use of the Girsanov theorem 
we obtain 



(64) 



By rearranging terms, the expectation on the right above can be rewritten in the form 



e *6 e E 'Io^-h E fIo^ ds 



e ±x 2 t+xEi Jo o- s ds jgC 



J*(Ei<T a +x)d£t-± Jo{E,u s +x) 2 ds 



(65) 



Since {^} is a Q-Brownian motion, it follows that the expectation appearing on the right 
side of (J55|) has the value unity. As a consequence, we deduce that 



(66) 



We see therefore that B t and H are independent, as claimed. We also see, by the form of 
its characteristic function, that B t is Gaussian, with mean and variance t. An argument 
similar to that presented above shows that if s < t then 



E P [ e xB s +y(B t -B s ^ =e ±^ e 



y 2 (t-s) 



(67) 



for all x, y, and hence Cov(B s ,B t ) = s for s < t. Since {B t } is a continuous Gaussian 
process and has the correct mean and autocovariance properties, we deduce that {B t } is a 
standard P-Brownian motion. 
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We have therefore shown that the solution to (0) can be put into the form (|52|). with 
£ t = H f*a s ds + B t , where H = lim^_ +00 (^ t |fl"|^ t ), and where {-£>*} is an independent 
standard Brownian motion. 

Conversely, as we have shown earlier in the paper, if on a fixed probability space (Q, J 7 , P) 
we are given a random variable H with the distribution ¥(H = Ej) = 7Tj together with an 
independent standard Brownian motion {B t }, then we can proceed as follows: we construct 
the process {&} by setting £ t = H f*a s ds + B t ; we define the process {H t } by H t = 
EfiJlJFj] where {Ft} is the filtration generated by {£ t }; we define the process {7^} by setting 
ix it = K[l{H=E i }\J~'t]] and we define the {J-tj-adapted Brownian motion {W t } by setting 

Wt = &- fo°sH s ds. Then {|^», defined by |^> = e^V*/ 2 ^ satisfies ©. 

Let us now consider briefly the "linearisation" technique (see, e.g., |8j and references cited 
therein) often used for studying the dynamics of ((TJ). Starting with (JTJ), we proceed as above 
to deduce defining the process {^} as in (jBT|) . Next we observe that if we introduce a 
process by defining 

\%) =exp(-iHt + ±H a s ^t~\H 2 cr 2 s ds^j ^0), (68) 
then \^ t ) satisfies 

d|tf t ) = -i#|tt t )dt - §a 2 # 2 |^ t >dt + \<r t H\%)ti£ t . (69) 

We note that \fy t ) appears in the numerator of (}52|) . Thus the relation between \ty t ) and 
|-0t) is given by 

IV*> = j^U , (70) 



where 



(%\%) = exp ^ (7 s i? s d6 - \ jf cr 2 # 2 ds) . (71) 

Now suppose we fix a finite time interval [0, T], and change to a new measure Q by use of 
the density $t = (^/tI^t) -1 - Then over the interval [0, T] the process {£t} is a Q-Brownian 
motion. Thus in the measure Q we can "solve" (pQ) by expressing |^) in terms of \^t), 
making use of (jHEJ) and f70|) . Equation (f6HJ) holds under both P and Q, but under Q it is a 
linear equation, and hence under Q we can regard ()68jl as the solution of ()69jl . 

The existence of this underlying "linearisation" of (0), together with the fact that there 
is only a single Hermitian operator appearing in the dynamics, may help to explain why the 
problem is exactly solvable. Of course, the "physics" is in the measure P, so any application 
of the linearisation technique to solve a physical model typically involves using the density 
Q^ 1 to change from Q back to P. In practical terms this means that realistic simulations of 
the trajectories of \ipt) cannot be efficiently achieved by use of the linearisation technique. 

On the other hand, our "random data" method involves a construction that is carried 
out entirely in the physical measure. The auxiliary measures are introduced simply for the 
purpose of verifying the results, not for the actual specification of the solution. In that sense, 
the random data method can be regarded as a major improvement over the linearisation 
method. In particular, to simulate a set of trajectories for {IV't)} we only need to simulate 
outcomes for H and {B t } in the measure P. 
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VII. VERIFICATION OF COLLAPSE PROPERTY 



In this section we shall verify directly that the solution ()48|) of the stochastic Schrodinger 
equation (JTJ gives rise to the collapse of the wave function. By substituting (jlj) into fT5|) . 
setting H = E k , and inserting the resulting expression into (f4*Kj) . we can express the solution 
of conditional on H = E k for some fixed value of A;, in the form 



J2i v^exp [~\Eit + \E t E k f* a 2 s ds + \E t f* a s dB s - \E 2 J* a s 2 ds) |&) 



TTi exp ( EiE k $1 a 2 s ds + E t f Q a s dB s - \E\ /* cx 2 ds 



1/2 



(72) 



In this situation we imagine that nature has "secretly" chosen the outcome H = E k (i.e. H 
takes this value for the given w £ O), and we want to show that the wave function evolves 
to the appropriate eigenstate. If we multiply the numerator and denominator of (J72j) by 
exp(— \E\ J * <J 2 s ds — \E k J^a s dB s ) and write uo ik = Ei — E k , then (|72|) becomes 



1/2 



It should be evident then, on account of condition (JHJ), that \ip t ) 
More precisely, defining 



1 V2 



■ c ) as t 



M t = exp f |a; ^ a s dB s 



aids 



.(73) 



oo. 



(74) 



we have, for any e > 0, 
P(M t > e) 



Flhw aABs-W cx 2 ds>lne 




> ^1 




(75) 



where iV(:r) is the standard normal distribution function 

1 



N(x) 



e 2 y2 d?/. 



(76) 
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Here we have used the fact that J* a s dB s /(J^ a 2 ds) 1 ^ 2 is normally distributed with mean 
zero and variance unity. We thus see that (JBJ) is satisfied if and only if P(M t > e) — > 
as t — > oo; and hence it follows that, given condition (jHJ), the state vector collapses to the 
designated eigenstate. The intuition behind this result is that the "signal" component of 
{rjt} eventually dominates over the "noise" component if © is satisfied. This is because the 
magnitude of J^a s dB s is on average about ( f Q ofds) 1 / 2 . 

We note, incidentally, that if the leading order behaviour of the integral of {of} is such 
that J *<j 2 ds ~ t a , with a > 0, then to leading order we have a t ~ t^ a ~ l \ and hence 
J * a s ds ~ t^ a+1 \ Since the magnitude of B t is on average of the order t 1//2 , we see that in 
this situation the signal component in (J3J also dominates over the noise component. 



VIII. REDUCTION WITHOUT COMPLETE COLLAPSE 

In this section we show that when condition © is not satisfied, state reduction neverthe- 
less takes place, in the sense that the energy variance decreases on average. However, unlike 
the models for which © is satisfied, in this case the terminal energy variance in general 
does not vanish. In other words, the state approaches an energy eigenstate, but does not 
get there. Physically, this situation corresponds to an approximate measurement of energy, 
in which some information concerning the energy of the system is revealed, but no definite 
outcome is obtained. 

To analyse this situation we consider the energy variance process flSBj) . which is given, 
equivalently, by 

V t = (^t\H 2 \^t) - (^t\H\^t) 2 . (77) 

Taking the stochastic differential of (J77|) and using the dynamical equation we find that 

dV t = -a 2 V 2 dt + a t K t dW t , (78) 

where K t is the third central moment of the energy: 

K t = {^ t \{H-H t f\^ t ). (79) 

We observe that the drift of {Vt} is strictly negative. Therefore, the energy variance is on 
average decreasing. However, if {at} is a square-integrable function, then {Vt} may converge 
to some finite nonzero value smaller than the initial value Vq. To investigate this scenario 
we write (J78|) in integral form: 

V t = V - [ a 2 K 2 ds + f o s kAW s . (80) 
Jo Jo 

Taking the expectation of each side, we obtain 

E\y t ] = V - [ a 2 s E[V s 2 ]ds. (81) 
Jo 

On account of Jensen's inequality we have E[V t 2 ] > (E[V t ]) 2 , and hence (JHTj) implies that 



E[V t ]<V - [ a 2 (E[K]) 2 ds. 
Jo 



(82) 
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Bearing in mind the fact that E[Vt] < E[Vs] for t > s, the inequality (J82j) implies 

E[V t ] < K,-(E[y t ]) 2 / a 2 s ds. (83) 
As a consequence, we obtain an upper bound on the expected value of the energy variance: 



nV t ] < ——j. — — \-l + Jl+Wo a*ds . (84) 




In the limit t — > oo the inequality ()84j) determines an upper bound for the asymptotic value 
of the expected energy variance. That is, on average the energy variance will be reduced to 
a value no greater than the asymptotic value of the right side of (|84|) . In particular, if (JOJ) is 
satisfied, then E[V t ] — > as t — > oo. On the other hand, if (jUJ) is not satisfied, then we can 
obtain a lower bound for E[Vy . Let V maiX denote the maximum possible variance that the 
energy can have, over all states. Then from (|81|) we get 

/•CO 

E[Ko] > V - KL / (85) 
Jo 

Hence providing J °° ofds < V /V^ ax we are ensured that state reduction will be incomplete. 



IX. FINITE-TIME COLLAPSE 

Having investigated the case in which the coupling {a t } decays too rapidly to lead to a 
complete collapse of the wave function, we turn to the situation where the integral of the 
coupling {a t } diverges over a finite time horizon. In particular, we consider the example 

* = (86) 

where a > and T > are fixed constants. For the resulting state vector dynamics we have 
the following stochastic differential equation: 



(rri \ 2 rri 
(H-H t ) 3 \rlH)dt+l^-^(H-H t )\rlH)dW t . 



(87) 



This model is of interest because the collapse of the wave function is achieved in finite time. 
The stochastic equation (|57j) is identical to the finite-time collapse model introduced in 
, where a solution to (}8*Tj) is obtained by considering an ansatz of the form 

C t =atH + (3 t . (88) 



11 



The noise term (3 t appearing here is a Brownian bridge t 28, 36] that vanishes at t — and 
at t — T. The vanishing of the noise at t — T guarantees the collapse of the wave function 
as t approaches T. In particular, since the coupling a in ()88|) is constant, {Q} is a Markov 
process, as is shown in |l2[. 
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In the present framework, it follows immediately from (jlj) that the appropriate ansatz 
for solving ()87|) is given by 

£ t = aT H In (J—^+B t , (89) 

where {B t }o<t<T is a standard Brownian motion. Remarkably, the two prescriptions (fKKj) 
and (|S9^1 give rise to the same solution to the stochastic equation (jHZJ). 

To see this we note that it follows from (j!6|) that the energy process associated with (JHTj) 
is given by 

E< ^ ex P (° TE * fa YZ-s ~ ¥ 2E i¥=l) 

Ht = 7^ yf-. (90) 

Ei *i exp [aTEi jj ^ d^ - 

In Section lllll we observed that under the probability measure Q the process is a 
Brownian motion. Therefore, in view of the expression appearing in the exponent of (|90jh 
we define a process according to the following scheme: 

Ct={T-t) f-i- d&. (91) 
Substituting (jSU|) into the right side of ()91|) we obtain 



(92) 



After a short calculation, we deduce that 



(T-t) f -J_d& = <rt# + A, (93) 
Jo 1 ~ s 



where {/3 t } is defined by 



fit = (T-t) I Y^- s dB s . (94) 



o 



However, we recognise in (J94j) a standard integral representation of a Brownian bridge [2 

jjj. It follows that {£"}, as defined by (JHIJ), 

can be put into the form ([88)1 . 

On the other hand, we also see that (}9~T]) is an integral representation of a Brownian 
bridge under Q, since in this measure {^} is a Brownian motion. Therefore, under Q, the 
energy process ()90|) can be expressed in terms of a Brownian bridge {£ t *} in the form 

„ = E,^exp(^ (aE^-^Eft)) 
1 E^ex P (^- t (aTE t a-^Eft))- [Jb) 

This result agrees with the result obtained in for the finite-time collapse model. 

We note, incidentally, that the innovation process associated with {£ t } in this example, 
given by 

W t = & - / a s H s ds } (96) 
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where a t = aT/(T — t), and the innovation process associated with {£ t *}, given by 

W t = Ct + J — ^ (C - <^TH^ ds, (97) 



as obtained in are identical if {/3t} is defined as in (I94j1 . This follows on account of the 
relation 



= £ + / Tf, Cds } (98) 

Jo T-s 

which can be verified by writing the right side of (|9*K|) in differential form: 
d£ + CAt = aHdt + dfc + (otH + fa) dt 

= eHdt - j^— t fodt + dB t + oHj^i dt + jr^- t Ptdt 

= H^- t dt + dB t = di t . (99) 

Therefore, the solution obtained here in terms of is equivalent to the solution obtained 



in 



111 ] using {£*}. The results above show that models exhibiting state- vector reduction over 



a finite time horizon are both feasible and tractable, and that such models can be usefully 
formulated by use of a time-dependent coupling. 
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APPENDIX A: PATH INTEGRAL APPROACH 

In this appendix we present a path integral method for calculating the conditional prob- 
ability process {7iu} and the energy expectation process {H t }. First we note that (J7|) can 
be written in the form: 

E[H\{Q < s < t ] = Y,EiTT it , (Al) 

i 

where n it = F(H = Ei\{^ s } < s < t ). By use of the Bayes theorem, the conditional probability 
is given by 



F(H = ^{Uck^) = ~„lT- tl iIT L - (A2) 





H = E i ) 


Si np({£.s}o< s <t 


H = E i ) 
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Here the expression 
p({£,s}o<s<t\H = Ei 



det(2vrS; 



exp 



J 



- Ei J " a s ds) (f B - Ei /J (x s ds) dwelt; )(A3) 



is the density function of the trajectory {£ s }o<s<t conditional on H taking the value E^. We 
write E(w, v) for the covariance of the random variables Ei f a s ds + B u and Ei a s ds + B v . 
The form of the density function (jA3|) follows from that fact that, conditional on H = Ei, the 
random variables {£s}o<s<t are jointly normally distributed. A straightforward calculation 
making use of well-known properties of Brownian motion shows that 

E(w,t>) = min(w,t>). (A4) 

To compute the inverse of the covariance we substitute (jA4|) into the relation 



/ E 1 {u, s)E(s, v )ds = 5(u — v ) 
Jo 



(A5) 



and differentiate the resulting expression twice in the variable v to obtain 

^- 1 (u,v) = -6"(u-v), (A6) 

where S"(t) denotes the second derivative of the Dirac delta function. 

We note that the conditional expectation of any functional of the trajectory {^ s }o< s << 
can be determined by use of the density function (|A3J) : the resulting expression corresponds 
to an infinite-dimensional Feynman integral. This is related to the fact that E(/u, v) gives 
the Feynman-Green function for a free particle [Til . 

To determine the form of the density function (|A3|) . one might consider substituting 



into (|A3|) and then applying integration by parts, since —8"(t) is a second-order differential 
operator (with the property that it has a positive "spike" at t = and a pair of negative 
"spikes" at t = + and t = 0~). However, the integral in the exponent of (jA3|) is ill-defined 
as it stands. Indeed, it is a straightforward exercise to verify that, depending on the order 
in which integration by parts is applied, one obtains different answers. To circumvent this 
difficulty we shall discretise the process {£ t } first, derive the corresponding expression for 
the conditional density function (jA2|) . and then take the continuum limit. 

Our strategy for determining the conditional probability law for the trajectory {^ s }o<s<t 
is thus as follows. Fixing t, we partition the range [0, t] into n equally-spaced intervals. 
In particular, we set Sk = kA (k = 0,1,..., n), where A = t/n. We then calculate the 
covariance of the random variables Ei J Q Sfe a u du + B. ,. and E^ J^ 1 a u du + B Sl , and obtain 
following expression for the covariance matrix (cf. |l6||): 
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The validity of this result should be evident from the continuous case (jA4|) . The inverse of 
the covariance matrix is the following second-order difference operator: 



Sl-SQ S2-S1 



V 



S2-S1 










1 s £ al l 







S2-S1 s 3 -s 2 s 3 -s 2 

J I + 








S3-S2 S3-S2 S4-S3 



S4-S3 



\ 






$n — 1 s n — 2 s n — l s n — 2 s n s n — l s n s n — l 

x — 



(A8) 
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Note that the last term on the diagonal is anomalous, and is different from the remaining 
terms on the diagonal. Because the partition of the interval [0, t] is equally spaced, we have 
s k — = A for all k = 1, 2, . . . , n. Therefore (jA8|) simplifies to 



yr 1 



A 



( 2 -1 0\ 

-12-10 

-1 2 -1 ■■■ 

... -1 2 -1 

0-12-1 

\ -1 1 / 



For simplicity of notation let us write 



oti(s k ) = £. Sk - Ei a u du. 

Jo 

Then the discretised form of the conditional density function (jA3|) takes the form 



(A9) 



(A10) 



p({£s 1 >U,---^s n }\H T = = ( — ) exp ( -lJ2J2 E ki l(X i( s khi( s i) 



k=l 1=1 



A v |n / . n-i 



2ty 



exp [ -£^2ai(s k )[a>i(s k+1 ) - ai(s k )] - ^af{s n ) J . (All) 



k=l 



Here we have used expression flA9|) for the inverse covariance matrix. Let us examine the 
terms in the exponent. Substituting definition (|A10|) we find that 



ai(s k )[ai(s k+1 ) - ai(s k )} = <U(£ Sfc+1 - <U) + / v u du / a v dv 

JO Js k 

(rsk+i fSk+i fSk 

(U+i-U) °udu-i Sk+1 l a u du + £ Sk a u du 

Js k J0 JO 



(A12) 



Turning to the conditional probability (jA2|) that we aim to determine, we observe by 
inspection of the right side of (jA2|) that all terms in the exponent that are independent 
of the eigenvalue E iy such as the term £s fe (£s fc+ i — £ Sfc ) in the right side of (jA12|) . cancel. 
This is because such terms appear in both the denominator and the numerator. Hence in 



23 



what follows we omit such terms. Equality modulo omitted terms will be denoted by the <■ 
symbol. Bearing this in mind, we see that the sum over k of the right side of (jA12|) gives 



n— 1 n— 1 



" x J* " r s k f s k+l 

) j ai(s k )[a i (s k+1 ) - ai{s k )) ~ EfjT o u du I a v dv 
k=i k=i Jo Js k 

+Ei ( XX(Ufx -&JA-& n j( o ^du + ^Aj . (A13) 

In deducing this result we have used the fact that for sufficiently large n, and hence suffi- 
ciently small A, the following relation holds to a high degree of accuracy: 

a u du = a Sk A. (A14) 
Similarly, omitting the term that contains no Ei we have 

a u d« + £ 2 f y a u dwj . (A15) 
Substituting (jA13|) and ()A15|) into (|A11|) . we deduce that 

(n-1 n-1 \ 

^^a Sfe fe fc+1 -6j-^ 2 ^<A . (A16) 
fc=0 fc=0 / 

Here we have made use of the fact that to a high degree of accuracy one has: 

/ ps n \ 2 n-1 n-1 

/ a u du) =^<(At) 2 + 2^a^A 2 . (A17) 

' k=0 k^l 

We have also used the fact that 

n— 1 fc— 1 n— 1 

k=l 1=0 k^l 

Note that the first sum in the exponent of (jA16j) begins from k — 0, and not k = 1 as in the 
right side of (|A13|) . This is because the last term in the right side of (|A13|) can be written 
as cr So 

(£si — £s )A (since ^ So = 0) and thus be absorbed in the sum. 
We are now in the position to take the limit n — * oo. In particular, the first term in the 
exponent of (jA16|) converges in this limit to an Ito integral of {at} with respect to the process 
{£t}, since the discrete approximation is always taken to be the value of the integrand at s k 
in each interval [s k , s k+ \]. The second sum, on the other hand, converges to the Riemann 
integral of the function {cr 2 }. Writing {pu} for the unnormalised density function given by 
the right side of (jA16j) . we thus deduce, in the limit n — > oo, that 

p it = exp (e { [ cr s d£ s - \E\ [ a 2 s ds) . (A19) 



As a consequence, for the conditional probability process {71"^} we obtain (fI5j). as desired. 
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APPENDIX B: DECOMPOSITION OF PATH INTO INCREMENTS 

In this appendix we introduce another method for obtaining {71"^}, based on the decom- 
position of the path of the quantum information process into its increments. The argument 
goes as follows. Recalling the dynamical equation (jHJ) satisfied by the information process 
{£t}, we note that {a t } moderates the strength of the signal showing the true value of H. 
An equivalent set of information can be obtained by moderating the noise level by use of 
the reciprocal function {l/a t }. Thus we consider a process {( t } defined by 

(t = Ht+ [ -dB s . (Bl) 

Jo °s 

The relation between {£ t } and {( t } is given by 

dCt = -d&. (B2) 



This alternative representation is motivated by Wonham |35(, where the filtering equation 
for a signal associated with a fixed random variable is investigated. 

The fact that the information implicit in {£ s }o<s<t and {Cs}o<s<i is equivalent will be 
shown. We begin with the discretisation of the period [0, t] into n equally-spaced intervals 
of size A = t/n, and set Sk = kA for k — 0, 1, . . . , n. We define a process for the increments 
of {Q} by setting y Sk = ( Sk+1 - ( Sk . It follows from (JBT|) that 

y Sk = HA + j o^B u . (B3) 

For each value of n, the random variables y Sk — HA (k = 0, 1, . . . , n — 1) are independent 
and normally distributed, with mean zero and variance 

a- 2 du. (B4) 

-h 

If we assume that A is small, then to a high degree of accuracy we have v Sk = cr~ k 2 A. Next 
we note that 

¥(H = E^, d, • • • , CJ = P (H = E^, y si , ■ ■ ■ , y^) . (B5) 

This relation follows from the fact that conditioning with respect to {Cs h }k=o,i,...,n is equiv- 
alent to conditioning with respect to the corresponding increments. Letting {vr 4 -"' ) } denote 
the conditional probability defined by (jB5|) . we conclude that 



7T 



exp (-ff:v-Hy Sk ~E t A) 2 ) 



it ~ / n-1 \ 

Ei n ex P -\ E v7* (Vs k - E i A ) 2 

V fc=o / 

(n— 1 n— 1 
EiA E v£y.„ - \E 2 £ v£X 
k=0 k=0 

/ n-1 n-1 \ 

E< *i E t A £ v-iy Sk - \Ej E ^ A 2 
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(B6) 
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Substituting the right side of (jB4|) into (|B6|) . and taking the limit as n gets large, we have 



In taking the limit to obtain a stochastic integral we have followed a line of argument similar 
to that of Appendix[X] This verifies that the information content of the trajectories {^ s }o<s<t 
and {C}o<s<t is the same. 
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